H2D+: a light on baryonic dark matter? 
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ABSTRACT 

' It has been suggested that the dark halos of galaxies are constituted by 
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cloudlets of cold (< 10 K) H2 and dense (> 10^ cm ^) molecular gas. Such 
gas is extremely difficult to detect, because the classical tracers of molecular gas, 



D ' CO and/or dust grains, have very low abundances and their emission is exceed- 

ingly weak. For this reason, the cloudlet hypothesis remains so far substantially 
^ ■ unproven. In this Letter we propose a new method to probe the presence of cold 

H2 clouds in galactic halos: the ground transition of ortho-Il2D+ at 372 GHz. 
> ■ We discuss why the 1120+ is abundant under the physical conditions appropriate 

for the cloudlets, and present a chemical model that predicts the H2D+ abun- 



•/^ ' dance as function of four key parameters: gas density and metallicity, cosmic 

■ ray ionization rate and dust grain size. We conclude that current ground-based 

^ ■ instruments might detect the ortho-H2D+ line emitted by the cloudlets halo, and 

' prove, therefore, the existence of large quantities of dark baryonic matter around 

CL|' galaxies. 



Subject headings: Cosmology - Dark matter 



■ 1. Introduction 

One of the most astounding results of modern cosmology is the fact that we currently 
firmly understand the nature of but a tiny fraction of the matter content of the universe: 
the "luminous matter" constituting the stars and gas in galaxies. This component amounts 
to only ~0.3% of the mass (Fukugita & Peebles 2004). The largest fraction of cosmic matter 
belongs to the so-called "dark sector" (95.4%), split into dark energy (72%) and dark matter 
(23%). In total, ~4.5% of the universe matter is present in the form of baryonic matter. 
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However, stars account for only ~5% of it. A significant fraction of the remaining 95% is 
thought to constitute the warm/hot intergalactic gas traced by the Lya forest (Savage et ah 
2002; Penton et ah 2004; Richtcr ct al. 2006) and/or X-rays (Cen & Ostriker 1999; Nicastro 
et al. 2005; Williams et al. 2006). However, the numbers are uncertain and up to 50% of the 
baryons at low redshift may still be unaccounted for. 

In galaxies, the evidence that most of the mass is not contained in the stellar or any 
other observed component stems from the study of rotation curves, derived from stellar light 
and HI gas emission (Rubin et al. 1962; Prugniel et al. 1998), which indicate the presence 
of a large invisible mass surrounding the galaxies. The nature of these dark matter halos 
has long baffled astronomers. It cannot be diffuse baryonic matter, because both neutral 
or molecular gas would be detected in one way or another. For this reason, it has been 
proposed that it may be non-baryonic in nature. However, the amount of dark mass in 
galaxies could be as large as the amount needed in the form of baryons in order to fulfill the 
constraints from Big Bang nucleosysthesis (Kirkman et al. 2003). Besides, there are other 
reasons to assume that the dark matter in galaxies may be mostly baryonic (Pfenniger et al. 
1994; Pfenniger & Combes 1994; Gerhard & Silk 1996). An invisible baryonic form would 
be condensed objects, like brown dwarfs or massive compact halo objects (MACHOS), but 
they both have now been ruled out as significant contributors (Alcock et al. 2001). Another 
possibihty is that baryonic dark matter is present in the form of cold (< lOK), molecular 
clouds, also called "cloudlets" (Pfenniger et al. 1994; Pfenniger & Combes 1994). There 
is considerable discussion in the literature about the allowed range of parameters for such 
clouds, and the values obtained depend on the particular model and physical constraints 
used. Gerhard & Silk (1996) consider both the collisional cross section of clouds and the 
thermal stability against collapse and conclude that the clouds should be colder than 10 K, 
have masses around 1 Mq and column densities above 1 x 10^^ cm^^. Walker & Wardle 
(1998) suggest that such clouds can be responsible for radio scintillations observed toward 
quasars and derive that the clouds should be a few AU in size and have masses below 
10~^ Mq. Wardle & Walker (1999) consider the thermal stability of such clouds against 
evaporation and find a lower limit for the mass of ~ 10~^ Mq to keep the clouds stable 
in the current epoch (low redshift). All these estimates involve uncertain parameters like 
the cosmic ray ionization rate in the halo. Rafikov & Draine (2001) use constrains from 
micro-lensing studies to derive a lower limit of the cloud masses of ~ lO'^Mp. and conclude 
that models with the Walker and Wardle parameters of 10 Mq and about 10 AU radii 
are still acceptable within the current observational limits. Kalberla et al. (1999) and Ohishi 
et al. (2004) suggest small halo cloudlets of 6AU and masses of 10~^ Mq as the explanation 
for the unidentified EGRET sources. It is, however, not clear if these clouds exist in a stable 
form, or if they are simply the highest density component of a fractal distribution of cold 
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molecular gas. In summary, several authors have invoked the possibihty that a substantial 
amount of baryonic mass is in cold cloudlets, but the uncertainty in the involved physics 
is so large that so far no consensus exists on the mass and size of the presumed cloudlets. 
Heithausen (2004) tried to probe the cloudlets in our own galaxy by means of high angular 
resolution CO observations. He indeed detected clumps of molecular material at 7-18 K, but 
it is unclear what their masses, sizes and densities are, given the uncertainty on the distance 
and CO abundance (see §2). In addition, Lawrence (2001) discusses the possibility that some 
faint SCUBA sources are in fact local very low temperature dust clouds, with characteristics 
similar to the halo cloudlets. 

In this Letter, we assume that the galactic dark lialos are constituted of such cloudlets, 
without further discussing the detailed structure and stability of such objects. Instead, we 
propose a new method to test this hypothesis: observations of the ground state transition 
of the ortho-H2D+. 

2. Why H2D+ 

Low metallicity cold molecular gas is intrinsically difficult to detect. The main con- 
stituent of cold gas, H2, does not have a permanent dipole moment and its ground transition 
can only be excited in warm (> lOOK) gas. The most sensitive probes of "standard" molecu- 
lar clouds in the galactic disk are dust continuum emission and CO millimeter line emission. 
However, if the gas is primordial (i.e. not enriched with the products of stellar nucleosyn- 
thesis) or very low in metallicity, the dust abundance is zero or very low and its emission 
undetectable. The same argument holds for CO. Even if the metallicity is non-zero and a 
small amount of CO and dust were present, the CO molecules would quickly freeze-out onto 
dust grains and in this way still escape detection. Similar chemical conditions can actually 
be found in particularly cold objects of our Galaxy: prestellar cores and protoplanetary 
disks, which contain regions so cold and dense, that CO and other heavy-element bearing 
molecules freeze out onto the cold grains. While galactic matter is still traceable through 
dust emission, the severe depletion of CO from the gas phase causes these regions to become 
"invisible" in terms of molecular emission, just like the presumed cloudlets in the dark galac- 
tic halos. However, recent breakthroughs in the observation and theoretical understanding 
of molecular deuteration have opened new ways to probe cold and CO depicted regions: the 
ground transition of H2D+ at 372 GHz, which is indeed a very specific signature of extremely 
cold and dense regions. As a matter of fact, the detection of this line has been sought after 
for decades in galactic giant molecular clouds (whose mass is > 10^ M©) without success 
(van Dishoeck et al. 1992), whereas it is easily detected in the small (< 1 M©) but cold and 
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dense pre-stellar cores (Caselli et al. 2003) and/or proto-planetary disks (Ceccarelli et al. 
2004). The initially unexpected presence of deuterated molecules in large abundances (Cec- 
carelli 2004) has now been explained by a new class of chemical models (Roberts et al. 2003). 
It is precisely the disappearance of the CO from the gas phase that is causing a dramatic 
enhancement of H2D+, HD^ and with respect to (for details see § 3). Excitingly, 
H2D+ and HD^ do have ground transitions observable with ground-based sub-millimeter 
telescopes (Stark et al. 1999; Vastel et al. 2004). These transitions can be used to probe the 
presence of the hypothetical cloudlets forming the baryonic dark matter. 



3. Model Description 

The chemical model wc arc using is largely identical with the one described by Ceccarelli 
& Dominik (2005), which is based upon work by Roberts et al. (2003) and Walmsley et al. 
(2004). We summarize here the basic assumptions entering into this model, and the modifi- 
cations we added to meet the requirement of cold gas clouds in galactic halos. In cold and 
dense gas, the most abundant molecule is always H2. HD is present as well, and generally is 
the main reservoir of deuterium, with an abundance of 3 x 10~^ relative to H2 (Linsky 2003). 
At low temperature, chemistry is generally driven by ion-neutral reactions, and therefore, 
by the process of ionization due to cosmic rays. Cosmic rays ionize H2 and quickly lead 
to the formation of H^, which is the main charge carrier in the gas. H^^ is destroyed by 
reactions with grains, heavy element bearing molecules like CO and N2 and, at a very low 
rate, through direct recombination with electrons. It also reacts with HD, and this is the 
starting point of deuterium enrichment chemistry whose overwhelming importance has only 
recently been recognized. H3 reacts with HD to form H2D+. 

H+ + HD ^ H2D+ + H2 (1) 

The reaction is endothermic, i.e. the reverse reaction H2D+ -|- H2 — > H^ -|- HD is practically 
forbidden at temperatures below ~ 30 K. Therefore, a deuterium atom captured in this way 
into an H3 isotopologue remains there until it is returned by a reaction neutralizing the 
ion. The relevant reactions are those with dust grains, CO and N2 molecules, and direct 
recombinations with electrons. As in cold and/or metal-poor environments the abundance 
of those species is low, H2D+ tends to include the entire positive charge. Of course, H2D"'" 
can itself react with HD to form HD^, and a further step leads to D^. 

Therefore, in cold gas, the positive charge is distributed between Hjj', H2D+, HD^, and D^ 
(Roberts et al. 2003; Walmsley et al. 2004). It turns out that, at extremely low metal- 
licities, this no longer holds. Under "normal" galactic conditions, H+ is quickly destroyed 



- 5 - 



by recombination on dust grain surfaces. At the lower grain abundances characteristic of 
low metallicity conditions, the only destruction channel for H"*" is direct recombination with 
electrons, a very inefficient reaction. Therefore, H"*" can dominate the positive charge. In the 
case where H+ dominates, the ionization balance can be calculated ignoring the presence of 
Hg and solving the equation 

xli+nn^kpe + 2[D]xH+nH2A;pHD - ^cr = (2) 

where Xfj+ is the H+ abundance with respect to H2, rin^ is the H2 density, [D] is the elemental 
deuterium abundance, .^cris the cosmic ray ionization rate, k^e and kpHD are the reaction 
rates of H+ recombination with electrons (3.6 x 10-'^'^{Tg^/300K)-^/'^ cm^ s'^; Millar et al. 
(1997)) and HD (2.0 x 10-%x.p[-A64/Tg^] cm^ s"^) respectively. In the case where H+ and 
its deuterated isotopologues dominate the positive charge, the electron abundance is found 
by solving Eq. 13 of CeccareUi Sz Dominik (2005). The low metallicity in primordial gas 
also affects directly the abundances of CO, N2 and dust grains. In the present model, the 
abundance of the grains is scaled with a factor Z/Zsun denoting the metallicity of the 
material relative to solar abundances. In the same way, we scale the abundances of CO and 
N2 linearly with the metallicity. 



4. Discussion and Conclusions 



Figure 1 shows the chemical composition of the cloudlets as function of four key parame- 
ters: the density in the cloud riu-i, ^^e metallicity Z, the cosmic ray ionization rate applicable 
in the halo, C,cr, and the size of dust grains present in the clouds, flgr. The plots show that 
H2D+ is almost never the main positive charge carrier in the cloudlets, even though it is 
almost always more abundant than Hj. The conditions are so extreme in the cloudlets that 
either D3 or H+ are the most abundant ions. As these are not observable, H2D+ remains the 
best tracer. It is very instructive to study the dependence of the H2D+ abundance on the dif- 
ferent parameters. First of all, the abundance deacreases as the density increases, leading to 
a weak dependence of the total H2D+ column density in a cloud on the cloudlet density. This 
also means that the internal structure of the cloudlets only weakly influences the resulting 
column densities. Also the dependence on the metallicity of the clouds is shallow. Changing 
it from 10^^ times solar to solar increases the abundance of H2D+ by just a factor of 30. The 
abundance of all ions scales approximately with cosmic ionization rate Finally, the 

size of grains plays an important role. The H2D+ abundance is highest for the smallest grain 
sizes. As the grains get bigger, the H2D+ abundance decreases but levels out at grain sizes 
of 0.1 /im. It should be kept in mind that the dependences discussed here are all relative 
to the chosen set of parameters. In reality, the parameter space is four dimensional and the 
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dependencies are complex. However, we can derive a range of abundances to be expected in 
the cloudlets. The H2D+ abundance is between 10~^ and 10~^^, when taking into account a 
large range for the four parameters varied simultaneously (larger than what is shown in Fig. 

The emission from the ground state transition of ortho-H2D+ at 372 GHz can be used 

to obscrvationally probe cold halo cloudlets. Like any other emission line, the ortho-H2D"'" 
line intensity depends on the excitation conditions (gas temperature and density), and on 
the average column density of ortho-H2D+ in the telescope beam. Since the critical density 
of the ortho-H2D+ ground transition is ~ 10^ cm^^, the line is very likely always thermally 
populated. In addition, the cloudlets temperature is likely lower than 10 K but larger than 
~ 5 K (the condensation temperature for H2), so that the intensity very weakly (less than 
a factor 4) depends on the temperature. Therefore, the most important parameter in the 
prediction of the H2D'^ line intensity is the ortho-H2D+ column density, which depends on 
the ortho-to-para ratio and the total H2D+ column density. 

a) H^D'^ ortho-to-para ratio: The chemical model gives the total abundance of H2D''". 
Theoretical estimates predict that for large (> 10^ cm"^) densities this ratio depends only 
on the gas temperature and reaches the unity at 9 K Flower et al. (2004). However, these 
theoretical estimates strongly depend on the ortho-to-para ratio of H2, itself a very poorly 
known (=observed) parameter. For the purpose of estimating the dctcctability of the ortho- 
H2D+ line at 372 GHz, we will consider an ortho-to-para ratio about unity. 

h) H^D^ column density: We assume that the halo contains a mass Mh in a sphere 
with radius i?h- The average H2 column density through this halo can be estimated to be 

^H2,h = Mh/nR'^fimp where = 2.8 is the mean molecular weight in H2 gas and rup is the 
mass of a proton. Note that this column is an average of lines of sight passing through a 
cloudlet, and other lines of sight not passing through one. In an uniformly distributed halo, 
the average column density of H2D'^ is then simply given by 

-^H2D+,h = Ti^rn ^H2D+ {nn2 , Z, ^cr, Ogr) (3) 

where ^HaD+is the abundance of the H2D+ molecule relative to H2. We will use the example 
of NGC3198 to have a guideline for the choice of general halo properties. This galaxy is 
considered the prototype for studies of halo dark matter. Interpretation of the HI rotational 
curve indicates that the halo extends more than 30 kpc, and the mass contained within 30 
kpc is 1.5 X 10^^ M0 (van Albada et al. 1985). If this mass is present as cold H2 clouds, 
it corresponds to an average total column A^Ha.h = 3 x lO^^cm"^. Using the results of the 
chemical modeling above, we arrive at an estimated column density of H2D+ between 3 x 10^ 
and 3 x 10^^ cm^^. 
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c) ortho-H^D^ line emission: If the ortho-H2D+ line is optically thin for an individual 
cloudlet, the line emission is simply proportional to the column density given above. If 
the line is optically thick, the emission is proportional to the column density at which the 
line becomes optically thick multiplied by the telescope filling factor. To have an order of 
magnitude, the ortho-H2D"'' line becomes optically thick at an ortho-H2D+ column density of 
~ 10^^, if the linewidth is ~1 km/s. Therefore, the line emission depends on the properties 
of individual cloudlets, which are very different depending on the particular model used (see 
§1). In the following, we will consider two cases. 

The first one assumes a cloudlet mass of IMq and a radius of 1000 AU (Gerhard & 
Silk 1996), corresponding to n^/j ~ 6 x 10'' cm~^ and an average cloud Nh^ ~ 1.3 x 10^^ 
cm~^. Using a line width of 2 km/s (twice the virial velocity of a single cloud), an abundance 
of x(H2D''")~ 3 X 10^^^ is the limit for optically thin emission. For the case of NGC3198, 
the fiUing factor of such clouds would be 4 x 10~^ when we assume a constant density of 
cloudlets throughout the halo. Thus, the maximum average ortho-H2D"'" column density 
of the halo is ~ 10^^ cm^^. For the second case, the least favorable to a H2D"'" detection, 
we use the unidentified EGRET sources which should have masses of 10~^ Mq and radii of 
10 AU. In this case, the cloudlet reaches 6 x 10^*^ cm~^ and the of such a cloudlet 
would be 2 X 10^^ cm~^. The ortho-H2D"'" line would be optically thick within the cloudlet for 
x(H2D"'") > 1.5 X 10~^^ (assuming that the linewidth is dominated by a 0.5 km/s turbulence). 
For an uniform dark halo mass and size hke in NGC3198, this case corresponds to a filling 
factor of 4 X 10~^, and, therefore, a maximum average ortho-H2D+ column density of ~ 10^° 
cm~^. 

Note that a centrally condensed distribution of the cloudlets could increase the filling 
factor and/or the average H2D+ column density close to the galaxy center significantly in 
both cases. For example, if the cloudlets follow a r"^ or power law the average column 
density in a 20" beam would be larger by a factor 4 and 60 respectively. In summary, the 
beam-averaged ortho-H2D"'" column density in the cloudlets halo of NGC3198, if uniform, 
is predicted to range between ~ 10^ and ~ 10^^ cm~^. If the halo is centrally condensed 
this number increases. For some cloudlet parameters used in the literature the emission will 
be optically thin, for others it will be optically thick. The two cases discussed above would 
suggest a ortho-H2D''" column density around to ~ 10^^ cm~^ or less, if the halo is distributed 
uniformly. If a power law applies, the ortho-H2D+ column density could be as high as 
~ 4 X 10^2 cm-2. 

At present, the rest frequency of ortho-H2D"'' at 372 GHz is observable with the CSO and 
APEX telescopes, and, for more distant galaxies, with JCMT, KOSMA and SMA. The fine 
detectability depends somewhat on the telescope employed, but, based on previous published 
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observations, it is safe to say that present facilities allow detections of ortlio-H2D+ column 
densities of about 1 x 10^^ cm~^ (CascIIi ct al. 2003; Ceccarelli et al. 2004). Much deeper 
integrations would lower this limit by about a factor 3, bringing the detectable ortho-H2D"'" 
column density to 3 x 10^^ cm~^. Therefore, we predict that the ortho-H2D"'" 372 GHz line 
should be detectable if the high end estimate of the H2D+ abundance in cloudlets applies. 
Indeed, a positive detection of this line would not only definitively prove the presence of 
important amounts of cold molecular material in galaxies, but would also strongly constrain 
the location and nature of the cloudlets. 

We warmly thank Francoise Combes for her frank and helpful discussions, and Dr. 
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We acknowledge Travel support through the Dutch/French van Gogh program, project VGP 
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Fig. 1. — Predicted abundances of the main charge carriers in the dark halo cloudlets (based 
on the Ceccarelli & Dominik (2005) model). Form left to right, the abundances are shown 
as a function of the gas density, metallicity, cosmic ray ionization rate and grain radius. 
When varying the metallicity, we reduced the standard abundances of dust grains, CO and 
N2 proportionally to the metallicity, relative to values typical in galactic molecular clouds: 
x(CO) = 9.5 X 10"^, a;(N2) = 2 x 10~^, and a dust-to-gas mass-ratio of 0.01. The plots have 
been done for a temperature of 8 K, density 10^ cm~^, metallicity 0.1, cosmic ray ionization 
rate 10""^^ s~^, and grain radius of 0.01 /im when these parameters are not varied in the 
relevant panel. The age of the clouds has been assumed larger than ~ 10'^ yr, which is the 
timescalc for the CO and N2 condensation onto the grains for a density of 10^ cm~^ (the 
timescale scales with the inverse of the density) . 



